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ABSTRACT 

Non-Gaussian velocity distribution in star forming region is repro- 
duced by inelastic clump collision model. We numerically calculated the 
evolution of inelastic hard spheres in sheared flow, which corresponds 
to cloud clumps in differential galactic rotation. This system fluctuates 
largely around equilibrium state, creating clusters with inelastic colli- 
sions and destroying them with shear motion. The fluctuation makes 
spheres have non-Gaussian velocity distribution with nearly exponential 
tail. How far from Gaussian distribution depends upon coefficient of 
restitution, which can produce the variety of degree of deviation from 
Gaussian among regions. 

Subject headings: Interstellar: Clouds - line profiles 
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1. Introduction 

In general, the velocity distribution in molec- 
ular clouds has an excess high velocity tail 
compared with Gaussian (e.g., b'algarone and 
IFhillips 1990t [Miesch and Scalo 1995|) . In par- 
ticular, Miesch and Scalo pointed out that 
the most important results are variety be- 
tween subregions of molecular cloud (in some 
subregions the distribution is nearly Gaussian 
while it deviates significantly from Gaussian 
in other subregions) and nearly exponential 
tail in the non-Gaussian cases ( [Miesch and] 
jcalo 1995[ ). 

Several explanations for the non- Gaussian 
distribution have been proposed: the high- 
velocity component came from unbound in- 
terclump medium ( [Blitz and Stark 1988| ) , 



hi ^h-speed clump collisions ( Ife to and fc grb 



tanzio i989[), nonlmear Alfven waves ([ET 



19901) , or the intermittency of interstellar tur- 



bulence ( [balgarone and Phillips 1990[ ). How- 
ever, it is difficult for these mechanisms to 
explain the variety and the exponential tail 
(Miesch and Scalo 1995). 

We propose an alternative mechanism, clus- 
tering of clumps through mutual inelastic 
collisions. This mechanism causes the non- 
Gaussian distribution with the variety and 
the nearly exponential tail with simple physics. 
Molecular clouds have hierarchical structure. 
Here we generically call density-contrasted ar- 
eas with any sizes "clumps". Clumps in- 
teract with one another through direct colli- 
sions and gravitational scatterings. In gen- 
eral, interparticle interactions in the exter- 
nal gravitational potential converts orbital 
energy to random motions (e.g., [Goldreich| 
and 'Iremaine 1982[ , [Ida 1990[ , phtsuki 1993[ ). 
Sheared motion between particles, which is 
due to differential rotation of galactic disk, 



causes encounters between them, even if they 
have no random motion. Thus random mo- 
tions are produced by interclump collisions 
and/or gravitational scatterings. On the other 
hand, interclump collisions also have effect to 
dissipate random motions, since they are su- 
personic, i.e., inelastic. Hence, some equilib- 
rium state of random motions of clumps are 
maintained. 

To make physics clear, we consider a simple 
model, inelastic hard spheres in sheared mo- 
tion and numerically calculated the dynam- 
ical evolution of the system. We neglected 
gravitational scattering, because what makes 
non-Gaussian distribution is inelaticity of col- 
lisions (see below) and the effects of gravita- 
tional scatterings are similar to those of elas- 
tic collisions (e.g., [Cuzzi et al. 1979| , phtsuki 

As we will show, the fluctuation around 
the equilibrium state is large in the case with 
relatively large inelaticity and densely packed 
spheres. That is, the mean value of the ve- 
locity dispertion fluctuates considerably and 
nonsteady spatial inhomogeneity of the dis- 
tribution of the spheres ( "cluster" ) is formed. 
This behavior is characteristic of a system 
with both energy input and dissipation and 
does not depend on details of dissipation 
mechanism. Hence, the essential features of 
realistic clump collisions can be described by 
our simple model, though dissipation is much 
more complicated in realistic clump collisions. 

The clustering of clumps produces non- 
Gaussian velocity distribution with nearly ex- 
ponential tail. The deviation from Gaussian 
is closely related to the degree of clustering, 
which is regulated by degree of inelasticity 
(restitution coefficients). In this model, the 
variety is naturally explained as different "ef- 
fective" restitution coefficients among subre- 
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gions. 

2. Method of numerical calculation 

Numerical setup of our model is as follows. 
We adopt rotating local Cartesian coordinates 
in a galactic disk. Hard spheres with normal 
restitution coefficient e are packed into two 
dimensional box having width Ly and height 
Lx- The X and y directions correspond to ra- 
dial and azimuthal directions of the galactic 
disk, respectively. We neglected "softness" 
of clumps, tangential restitution coefficient, 
and the thickness of the disk for simplicity, 
which hardly affects the properties of the re- 
sults here. 

For the limitation of calculation, number of 
spheres, i.e., and Ly, cannot be taken so 
large. Hence the box corresponds to the area 
which is much smaller than the area identified 
to one subregion by observation. We describe 
dynamics of a whole subregion by time aver- 
age of this box instead of by summation of 
such boxes which cover a whole subregion. 

This box has periodic boundary condition 
in the x direction and sheared boundary con- 
dition in the y direction, corresponding to 
sheared motion in a galactic potential (Fig.|I|). 
When a sphere passes through lower horizon- 
tal side of the box downward, it re-enters 
into the box from upper horizontal boundary. 
However, horizontal position differs between 
exiting position and entering position. Sup- 
pose that y coordinate of the position where 
the sphere goes out of box is yout and that of 
entering is yin. Then is equal to yout + Ut, 
where and t is time and U is sheared velocity 
between positions at upper and lower bound- 
aries. If yout + Ut is outside the box, y^ is 
taken equal to be yout + Ut — Ly. The y com- 
ponent of velocity after re-entering box also 
increase by U;Vin = Vout + U where f j„ and 



Vout are the y component of velocity when 
the sphere exits from and re-enters the box 
respectively. When a sphere passes through 
horizontal boundary upward, yout should be 
yin — Ut and Vin = Vout — U . On the other 
hand, if a particle passes horizontally the ver- 
tical boundaries, Xin = Xout and Uj„ = Uout, 
where u is the x-component of velocity. 

In addition to self-gravity of clumps, we ne- 
glected Coriolis force and tidal force (the dif- 
ference between galactic potential force and 
centrifugal force) in the equations of motion, 
because mean collision time is much shorter 
than galactic rotation time, i.e., characteris- 
tic time of Coriolis and tidal forces. (Prelim- 
inary results by 3D A^-body calculation with 
self-gravity, Coriolis and tidal forces show es- 
sentially the same results.) 

3. Results and their physical interpre- 
tation 

By using this model, we perform numer- 
ical simulations. 256 spheres with diameter 
1.0 are considered. The box size L^. x Ly is 
taken to be 18.2 x 15.2 and sheared velocity 
U = is 0.5. The box size is taken to realize 
dense packing (in this case ~ 70%). Clusters 
of clumps are created by inelastic collisions 
and are smeared out by sheared motion. If 
packing is so sparse that mean collision time 
is longer than smearing timescale, clusters are 
not maintained. Hence dense packing is re- 
quired to produce non-Gaussian distribution 
due to clustering. The choices of values of L^, 
Ly, and U do not change results, as long as 
dense packing is realized. 

Restitution coefficient e is taken to be 0.6 
and 0.9. Result is sensitive to value of e. Af- 
ter discarding transient period starting from 
initial condition with equally spacing and ran- 
dom velocity (g [—0.05,0.05]), we measure 
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distribution of velocity component over 
120 (180) snapshots for e = 0.9 (e = 0.6) 
every 256 (= total number of spheres) colli- 
sions (Fig.0)0. The distribution in Figure§(a) 
clearly deviates from Gaussian distribution. 
It has excess high-velocity tail that is nearly 
exponential. This property is independent of 
initial conditions when e is relatively small. 
This universal property is consistent with ob- 
servations. The degree of the deviation in- 
creases with deviation of the restitution coef- 
ficient e from unity. The result with e = 0.9 
shows almost Gaussian distribution (See Fig. 
0(b)). 

We found that non-Gaussian distribution 
is made by unsteady clustering due to in- 
elastic collision. In Fig. ^, time evolutions 
of mean velocity dispersion are shown. The 
vertical and horizontal axes are mean veloc- 
ity dispersion and cumulative collision num- 
bers. The velocity dispersion appreciably 
varies with time in the e = 0.6 case. In the 
e = 0.9 case, the variation is smaller. In the 
limit of e — *• 1, the velocity dispersion fluc- 
tuates only within statistical one. We found 
that this variance is induced by clustering of 
clumps. Generally speaking, introduction of 
inelastic collision causes spatial inhomogene- 
ity. Goldhirsch and co-workers ( |GoldhirscIi 
an( |i Zanetti 1993| , [Goldhirsch et al. 1993|) 
have already pointed out that spatially ho- 
mogeneous distribution of hard spheres with 
inelastic collisions are unstable and collapse 
into several clusters. In our case, shear mo- 
tion destroys the clusters. Hence creation and 
destruction of clusters are repeated. Through 
these processes, energy input from shear mo- 



^This is because distribution of v component is 

strongly deformed by shear velocity 
^Thus total number of collisions observed is 120x256 — 

30720 for e = 0.9 and 46080 for e = 0.6. 



tion is intermittently dissipated by inelastic 
collisions, which makes large variation of ve- 
locity dispersion. 

This time variance of dispersion does make 
the deviation from Gaussian distribution, since 
a composition of many Gaussian-like distribu- 
tions with different dispersions is non-Gaussian 
distribution with excess high-velocity tail. Thus 
non-Gaussian velocity dispersion can be ex- 
plained by nonsteady clustering due to inelas- 
tic collisions. 

As noted before, the time variance of the 
present calculation corresponds to spatial vari- 
ance in subregion of molecular cloud. Hence, 
the observed non-Gaussian distribution with 
excess high-velocity tail can be explained by 
spatial variance of clustering within a subre- 
gion caused by inelastic interclump collisions. 

This model explains the observed variety 
of the distribution as the variety of "effective" 
restitution coefficient. As shown in Fig. 0, the 
deviation from Gaussian depends upon resti- 
tution coefficient. Effective restitution coeffi- 
cient would be different among subregions in 
clouds. For example, effective e can change 
with p as follows. The relative importance 
between inelastic collisions and gravitational 
scatterings is indicated by the ratio of sizes 
r of clumps to their tidal (or Roche-lobe) 
radii rt defined by {GM/^l'^y^^ where M is 
clump mass and Q is galactic rotational fre- 
quency; if r/rt <C 1, gravitational scatterings 
are more important, and vice versa (e.g., Oht- 
suki 1993). For clumps of molecular cloud in 
the neighborhood of the sun, we obtain 



^^.0.1 ^ 

n I1OM0PC-2 



^1/3 



where p is mass density of clump. This 
means that direct collisions effective as 

gravitational scatterings in evolution of ran- 
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dom motion of clumps, analogous to plane- Ohtsuki, K. 1993, Icarus, 106, 228. 

tary rings. Since gravitational scattering is 

equivalent to e = 1 collisions, effective e can 

change with p through Eq.(l). Furthermore, 

different configurations and random velocities 

of clumps may change effective e. 

Thus inelastic interclump collisions explain 
the nearly exponential high-velocity tail and 
the variety of distribution, and hence it can be 
a basic mechanism why velocity distribution 
in molecular cloud deviates from Gaussian. 

We would like to thank Hideki Takayasu 
for fruitful discussion. 
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Fig. 1. — Local Cartesian coordinates with 
sheared boundary condition. 

Fig. 2. — Velocity distribution (a) e = 0.6 
(b) e = 0.9 w is renormalized so as to have 
variance of unity. 

Fig. 3. — Time evolutions of mean velocity u 
dispersion are shown. Upper: e = 0.9, Lower: 
e = 0.6. Data before 5120th coUison are dis- 
carded. 
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